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The continuous quantum phase transition between noninteracting, time-reversal symmetric topo¬ 
logical and trivial insulators in three dimensions is described by the massless Dirac fermion. We 
address the stability of this quantum critical point against short range electronic interactions by 
using renormalization group analysis and mean held theory. For sufficiently weak interactions, we 
show that the nature of the direct transition remains unchanged. Beyond a critical strength of 
interactions we find that either (i) there is a direct first order transition between two time reversal 
symmetric insulators or (ii) the direct transition is eliminated by an intervening time reversal and 
inversion odd “axionic” insulator. We also demonstrate the existence of an interaction driven first 
order quantum phase transition between topological and trivial gapped states in lower dimensions. 
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Introduction : The spin-orbit coupled, time reversal 
symmetric insulators in two and three dimensions be¬ 
long to the Altland-Zirnbauer class All mm- Due to 
the existence of nontrivial Z 2 topological invariants in 
both spatial dimensions, a band insulator can be clas¬ 
sified either as a strong topological insulator (TI) or a 
trivial/normal insulator (NI). The four component mas¬ 
sive Dirac fermion provides an efficient low energy de¬ 
scription of Kramers degenerate conduction and valence 
bands in such systems EHS]- The sign of the Dirac mass 
provides the topological distinction between two insulat¬ 
ing states. In a clean, noninteracting system the univer¬ 
sality class of a continuous topological quantum phase 
transition (QPT) between them (where the Dirac mass 
vanishes) is described by a massless Dirac Hamiltonian. 
Recently, there has been considerable interest in under¬ 
standing the nature of such topological transition in the 
two-dimensional HgCdTe quantum well [7] and three di¬ 
mensional materials such as BiT^Si^Se^ (5] [2] and 
(Bii_ x In x ) 2 Se 3 [0. IT]. Although the phase diagram 
of noninteracting, disordered systems has been investi¬ 
gated in many analytical and numerical works |12l42()j 
and the interplay of disorder and long range Coulomb in¬ 
teractions has also been studied perturbatively for both 
two m and three dimensions [TB], the effects of strong 
short range electronic interactions on the phase diagram 
of clean All insulators is not well understood. The inves¬ 
tigation of this fundamental problem is the central theme 
of this Rapid Communication. 

We show that the three-dimensional massless Dirac 
fermion is stable against sufficiently weak, but generic 
short range electronic interactions. Consequently, the 
universality class of the continuous QPT between the 
TI and NI remains unaffected up to a critical strength 
of interactions, beyond which the quantum critical point 
(QCP) becomes destabilized due to the spontaneous gap 
generation for the massless Dirac fermion. There are 
two possibilities for the gap: (i) an inversion ( V ) and 
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FIG. 1: Zero temperature phase diagram of three-dimensional 
Dirac fermion for (a) g w = 0 and (b) g a = 0, for b = 0.2. All 
the parameters are dimensionless (defined in the text). g a 
and g, r are interaction couplings for scalar and pseudoscalar 
mass generation, respectively, and m is the band mass for 
Dirac fermions. Dirac semimetal is realized on the solid line. 
The transition between two insulators is first order (continu¬ 
ous) across the dotted (solid and dashed) line(s). The three 
phases TI, NI and AI meet at a multicritical point as shown in 
(b). For d=3, similar first order transition can occur between 
strong and weak topological insulators, and also between weak 
topological and normal insulators. 


time reversal (T) preserving scalar mass, which does 
not break any boncifide microscopic symmetry and (ii) 
a pseudoscalar mass that spontaneously breaks genuine 
discrete V and T symmetries and corresponds to an “ax¬ 
ionic insulator” (AI), which has recently been discussed 
in the context of magnetic insulators EDH2I, a Kondo 
singlet phase [23), and p + is superconductor [24] . The 
scalar mass generation leads to a direct first order QPT 
(without band gap closer) between the TI and NI [see 
Fig. l(a)| . This is a fluctuation driven first order transi¬ 
tion and is distinct from the generic first order transition 
between two different broken symmetry phases 1251 !f&f . 
On the other hand, the nucleation of pseudoscalar mass 
separates two T symmetric insulators, while eliminating 
a direct transition between them [see Fig.|l(b)|, and three 
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phases (TI, NI and AI) meet at a multi-critical point 
(MCP). In contrast to the TI and NI, which are respec¬ 
tively characterized by the quantized magnetoelectric co¬ 
efficients 7 r and 0, the AI phase supports a nonquantized 
magnetoelectric coefficient. Consequently, the AI phase 
does not possess gapless surface states, and the massive 
fluctuations of the magnetoelectric coefficient cause dy¬ 
namic magnetoelectric effects. We also reach similar con¬ 
clusions regarding the fate of the QPT between the TI 
(quantum spin Hall insulator) and NI in two dimensions. 
However, the nature of the possible T breaking phases in 
two and three dimensions are considerably different. We 
also show a similar first order transition between topo¬ 
logical and trivial superconductors in one dimension. 

Model: We will mostly be interested in the three- 
dimensional materials with V and T symmetries. For 
concreteness, we choose the following tight binding 
Hamiltonian on a cubic lattice jT 

3 

H = £*Lk£ Tj sin(fcja) + T 4 M]'P k 
k j =i 

3 

+t 2 X! £ ^k r 4 t 1 - COs(fejO)] 'Ffc, (1) 
k i=1 

as a minimal model for describing the insulating states 
in class All, where a is the lattice constant. The four 
component spinor Tj = (c +;t)k , c +4ik , c_ itjk , c_ )4 . ik ) is 
comprised of the electron annihilation operator c riS-k for 
states with wavevector k, the orbital parity r = ± and 
the spin projections s =t / i- The relevant discrete sym¬ 
metry operations are defined according to V: k —> — k 
and *F k —)■ r 4 >F_ k , and T: k — > — k and \F k —> TiTa'I'^. 
We have used only four mutually anticommuting matri¬ 
ces Tj = cj j g) ri (with j = 1,2,3), T 4 = cto < 8 > 73 in 
Eq. ([!]) to enforce simultaneous V and T symmetries. 
The bilinear ’F'Ts’F, where Fs = ao < 8 > t 2 , maintains two 
fold degeneracy of the valence and conduction bands, but 
separately breaks V and T. By contrast, fermion bilin- 
ears 'F'f£ a b'F, with £ a h = [T a ,Ph]/(2i), either break V 
or 7", and lift the Kramers degeneracy. Consequently, 
’F 1 '£ ab ’F and are absent in H. Here, a^ and 

are two sets of Pauli matrices respectively operating on 
the spin and the parity indices. The above model de¬ 
scribes both strong (for —2 < ^ < 0 ) and weak (for 
—4 < |^ < —2) TIs. The strong TI and NI phases re¬ 
spectively correspond to sgn(Mt 2 ) < Oandsgn(Mf 2 ) > 0 
and M = 0 describes the QPT between them. In con¬ 
trast, the transition between weak TI and strong TI (NI) 
takes place when ^ = —2(—4). For simplicity, we focus 
on the strong TI-NI QCP. But, our results are equally 
applicable near the weak Tl-strong TI or NI QCPs. 

Since the minimal gap in the spectrum near the strong 
TI-NI QPT occurs at the T point k = (0,0,0) [27], the 
low energy quasiparticles in its vicinity are described by 


the following continuum Hamiltonian 


H 3d 


d 3 k + 

- *F 

(2tt) 3 k 


3 

v ^Tjkj +r 4 (M + Bk 2 ) 

j=l 


'bk, 


(2) 

where v = t±a and B = t 2 a 2 . The integral over the mo¬ 
mentum is restricted to an ultraviolet cutoff A. The TI 
and NI phases respectively correspond to sgn (MB) < 0 
and sgn (MB) > 0. At the QCP (M = 0) between these 
two insulators the critical excitations are described by 
massless Dirac fermions and Bk 2 act as a momentum de¬ 
pendent Wilson mass. Only the fixed point Hamiltonian 
of massless Dirac fermion (M = B = 0) possesses a global 
chiral U(l) symmetry under which *F —> 4>e ieS45 [2S]. At 
the microscopic level B / 0 and H 3d enjoys only a re¬ 
duced Z 2 particle-hole symmetry, captured by the anti¬ 
commutation relation = 0. Additionally, af¬ 

ter an appropriate modification of the spinor, the same 
Hamiltonian operator can also describe the topological 
QPT (BEC-BCS transition) for superconducting systems 
in Altland-Zirnbauer classes AIII and DHL 


For simplicity, we here restrict ourselves to class All 
and consider the effects of repulsive, short-range inter¬ 
actions. For a V, T symmetric system with rotational 
symmetry, generic model of short-range interaction 


‘ l) r; f — 


er x 


y (^roT) 2 + ^ (^£45vF) 2 


+Y ('F t r 4 'F ) 2 + y ('F+rs’F ) 5 


( 3 ) 


is described by only four linearly independent couplings 
due to the Fierz identity [29]. For Ai = A 2 = 0 and 
M = B = 0, H 3 d + H int corresponds to the celebrated 
Nambu-Jona-Lasinio (NJL) model for mass generation 
through spontaneous U(l) chiral symmetry breaking [501 . 

RG analysis: Since the dynamic scaling exponent for 
massless Dirac fermion is z = 1, the scaling dimensions 
of the quartic interactions are [Aj] = {z — d) = (1 — d). 
Therefore, any sufficiently weak short range interaction 
is an irrelevant perturbation at the massless Dirac QCP 
for d > 1, and leaves the universality class of the direct 
TI-NI transition unchanged. In the weak coupling limit, 
AiS only modify the phase boundary in a non-universal 
manner. These features and the potential breakdown of 
the massless Dirac theory for strong interactions can be 
captured through a renormalization group (RG) analysis, 
controlled via simultaneous 1 /Nf- and e-expansion about 
the lower critical dimension d = 1, where Nf is the flavor 
number of four component fermions and e = d— 1. Within 
the Wilsonian momentum shell method, we integrate out 
the degrees of freedom inside —00 < oj < 00 and Ae _i < 
k < A, and subsequently rescale according to x —> xe l , 
t ->• re\ 'F e~ dl / 2 ^. In the N f —> 00 limit, we obtain 
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the following flow equations 


which takes the form 


dm 

hi 


m( 1 + 2 g a ) + 2 g a b, 


db 

dl 


dg 2 dg a 

dl ~ e92 ’ dl 


-eg a + 2 3 ", 


-6, 

dfU 

dl 


dgi 

dl 


= ~eg i, 


— —effir + 2g 2 , (4) 


where we have introduced the dimensionless couplings 
9i = Ai5' d 7V/A d_1 /[u(27r) d ], m = M/(v A), & = BA/v, 
and Sd is the surface area of a d-dimensional unit sphere. 

For sufficiently weak interactions, Dirac mass m is the 
only relevant variable. In this regime, gi(l) ~ gt, oe _eJ , 
b(l) ~ boe~ l , where gi : o and &o are the bare values of 
the corresponding couplings and the solution of Eq. Q 
yields 


i(0 + ( ^y L K05a(0 = e ; + ^-j-&oS<r,o^ 


( 5 ) 


Therefore, m* = Too + acts as an effective mass 

and the TI-NI phase boundary is determined by m* = 0. 
This agrees well with the phase boundary determined 
through minimizing the free energy for weak interactions, 
as shown in Fig. |f (a)| for d = 3 [3j. The flow equations 
also show that g\ and g 2 are always irrelevant pertur¬ 
bations, and can be ignored. Thus, we can restrict our¬ 
selves to the NJL model that supports three unstable 
fixed points: (i) g c G = e/2, g° = 0; (ii) g% = 0, g* = e/2; 
and (iii) gf = g% = e/2. The first two are QCPs, re¬ 
spectively describing the generation of scalar and pseu¬ 
doscalar masses for massless Dirac fermion at strong cou¬ 
plings, whereas the third one is a bicritical point. For 
d > I and at strong coupling (g.; > g /), the appropri¬ 
ate correlation length £, ~ A~ 1 (g i — g < [)~ IJi provides the 
infra-red cutoff for the RG flow, with a correlation length 
exponent Vi = 1/e. Therefore, for d — 2 and d = 3 we re¬ 
spectively have non-Gaussian and Gaussian (mean-field) 
itinerant QCPs [32J. As b(l) = boe~ l , we can ignore it in 
the RG sense in the weak-coupling regime. However, the 
effects of b cannot be neglected for gi > gf s. For eluci¬ 
dating the dramatic effects of b in determining the nature 
of the strong coupling phases and the associated transi¬ 
tions, next we adopt the method of large Nf mean-field 
theory. 

Free energy. We perform the Hubbard-Stratonovich 
decoupling of the quartic terms proportional to g a and 
g n , with bosonic fields E and n, respectively, that cou¬ 
ple to the fermion bilinears as E and n TYs'I'. 

When the corresponding coupling constants exceed their 
critical strengths, these bosonic fields acquire expectation 
values, while giving rise to the gap in the spectrum. With 
(n) = 0, (E) ^ 0, (E) + M acts as the effective scalar 
mass, with TI and NI being the only possible phases. In 
contrast, when (n) ^ 0, V and T are spontaneously bro¬ 
ken, giving rise to the AI phase. Since the free energy 
density (F) has the dimension EL~ d ~ L~^ d+1 \ we de¬ 
fine a dimensionless quantity / = F{2n) d /(NfSdvA d+1 ), 


2 2 /*1 

/ =-^—-2 / dxx d ~ 1 hx 1 + (cr + to + bx 2 ) 2 + n 2 . 

T J q 

( 6 ) 

We have also introduced the following dimensionless vari¬ 
ables a = E/(uA), 7 r = n/(uA), x = k/A. 

Note a = —m, n = 0 corresponds to the massless 
Dirac fermion describing the continuous QPT between 
the TI and NI, and remains as the global minimum of / 
up to critical strengths of ^s. Consequently, the phase 
boundary between the TI and NI is determined by 


m = 


2g* 
3 b 3 


(2 — 6 2 )\/l + b 2 



( 7 ) 


in d = 3, which is equivalent to m* = 0 for small b , up 
to the critical strengths of the quartic interactions 


9 a = 


1 - 
2 . 


1 + b 2 + y/l + b 2 


9l = \ 1 + Vl + b 2 


■ ( 8 ) 


Notice that g* a > g* for an arbitrary value of b. This is 
related to the fact {r^rs} = 0. Thus a finite b favors 
the nucleation of the pseudoscalar mass. In a continuum 
model without any momentum dependent mass (b = 0) 
the critical couplings g* = g* = e /2 = 1 become identical 
to the ones found from the RG analysis of the NJL model. 

We numerically minimize the free energy to obtain the 
phase diagrams in Figs. |l(a)] and |l(b)| However, all the 
salient features can be understood by expanding / for 
small cr + to and 7r, and retaining only the lowest order 
^-linear contributions. The free energy is then given by 


fzd 

' 1 
.9 7T 


m 2 1 


m b 

1 

cr 2 

' 1 1' 

hdcr ~ 2 

— (7 

.9a 2 

+ Y 

.9a 9%. 


b 2 2 \ (o' 2 + 7r 2 ) 2 

-cr(cr + 7T ) + A---log 


2e -1 / 4 


( 9 ) 


for d = 3, after shifting a —b cr — m. When we consider 
the nucleation of the scalar mass (for sufficiently strong 
g a ), the following subtleties need to be addressed: (i) 
(^ + |) acts as the external field for a and (ii) the 
momentum dependent mass oc b gives rise to all odd 
powers (cubic and higher) of cr in the expression for the 
free energy density. Due to the presence of all the odd 
powers of a, the nucleation of the scalar mass proceeds 
through a first order QPT between the TI and NI [see 
Fig. l(aff . By contrast, the pseudoscalar mass genera¬ 
tion (for sufficiently strong g v ) occurs through continu¬ 
ous QPT. Three phases (TI, NI and AI) meet at a MCP 
(located at m* = 0 and g n = g* in the g a = 0 plane), 
as shown in Fig. |l(b)| The MCP displays mean-field ex¬ 
ponents, but the hyperscaling is violated by logarithmic 
corrections [due to the 7r 4 log(7r 2 ) term in Eq. ([£])]. 

Two-dimensions: Our study can be generalized to ad¬ 
dress the transition between the TI (quantum spin-Hall 
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insulator) and NI in two dimensions. In the continuum 
limit, these two insulators and the transition between 
them can be described by the Bernevig-Hughes-Zhang 
model [53] 

H 2d = J t[v(r 3 k x +T 5 k y )+T 4 (M + Bk 2 )]tf k . 

( 10 ) 


The TI and NI phases in this model are realized for 
sgn (MB) < 0 and sgn (MB) > 0, and the continu¬ 
ous QPT between them (M = 0) is described by two- 
dimensional massless Dirac fermion. Strong interac¬ 
tions can give rise to the following four Dirac masses: 
E = ('P'iT 4 'P), n lj2 = (4' t ri 2 'h), and 5 = (^T, 35 ^). 
As in d = 3, the nucleation of the scalar mass E (for suf¬ 
ficiently strong g a ) proceeds through a first order QPT, 
which can be understood from the following free energy 
density 


f2d — 


rrc 2 
2 9a 3 

2 


. m 2 

~a[ — +~b 

9a 3 


a 

y 


+ M‘ o + b sgn(cr) + 0(b , a ) 


( 11 ) 


where g c a = 1/2. The corresponding phase diagram is 
qualitatively similar to Fig. 1(a) In contrast, conden¬ 


sation of T-odd magnetic masses rEx 2 or the anomalous 
charge Hall mass S [33] takes place through a continuous 
QPT as in Fig. |l(b)| (see also Ref. [55]), giving rise to 
a MCP in the (m, g v ) plane (with non-mean-held expo¬ 
nents [36] [37]). Our proposed first order transition be¬ 
tween TI and NI has recently been observed in numerical 
works in both d = 2 [35] and d = 3 [33]. 

One dimension: In one dimension, the quartic interac¬ 
tion (g a ) is marginally relevant and destabilizes the mass¬ 
less Dirac fixed point for infinitesimal strength [setting 
e = 0 in Eq. 0] [40] . In the presence of a momen¬ 
tum dependent mass, the QPT between topological and 
trivial gapped states can be first order for a sufficiently 
large number of flavors. This is the reason behind the 
existence of a fluctuation driven first order QPT for the 
TV-color Aslrkin-Teller chain when N > 3 [4TU431 . The 
pertinent continuum model is described by N species of 
two component Majorana/Jordan-Wigner fermions with 
a k dependent mass in the vicinity of the decoupled Ising 
QCP, and microscopic four-spin coupling gives rise to 
0(N ) invariant quartic interaction. Such a first order 
QPT can be germane for describing the direct QPT be¬ 
tween topological and trivial superconductors in different 
classes. 

Conclusions: Our results can also shed light onto the 
finite temperature phase diagram of strongly interacting 
Dirac fermion. At finite temperatures and for sufficiently 
strong interactions, the massless Dirac fermion can un¬ 
dergo either (i) a first-order classical phase transition and 
enter into the TI or NI phase, or (ii) a continuous classical 


phase transition, giving rise to a T-breaking insulator. In 
general, the semiconductors (where most of the TIs have 
been identified) are weakly correlated materials. How¬ 
ever sufficiently strong electronic interactions (local and 
nonretarded) can be mediated by optical phonons, below 
the scale of optical frequency [33115(11 . 

Finally, we comment on the effects of disorder on var¬ 
ious QPTs, discussed for clean systems. Harris criterion 
dictates that any continuous QPT is stable against weak 
disorder if v > 2/d [37] and this has been explicitly shown 
for the TI-NI QPT at d = 3 (T6] . Since u = 1/2 < 2/3 
at the clean MCP in d = 3, the universality class will be 
changed by an infinitesimal amount of randomness into 
a disorder controlled class satisfying v > 2/3 [35] , At 
present it is not well established whether the exact value 
of v for MCP in d = 2 is bigger or less than one (= 2/d), 
and we cannot properly assess its stability against weak 
disorder. On the other hand, in the thermodynamic 
limit, the first order transition (both classical and quan¬ 
tum) will be rounded by disorder into a continuous one 
in d = 1, 2, while it can survive in d = 3 for sufficiently 
weak randomness [391155] . The analysis of the universal¬ 
ity class at a putative disorder and interaction controlled 
critical point is beyond the scope of the present Rapid 
Communication. 
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